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Abstract 

We study the one-dimensional Dirac equation with local PT-symmetric 
potentials whose discrete eigenfunctions and continuum asymptotic eigen- 
fimctions are eigenfunctions of the VT operator, too: on these conditions 
the bound-state spectra are real and the potentials are reflectionless and 
conserve unitarity in the scattering process. Absence of reflection makes it 
meaningful to consider also PT-symmetric potentials that do not vanish 
asymptotically. 

1 Introduction 

Reflectionless potentials, i.e. potentials that are transparent to incident waves 
at all energies, have played a special role in quantum mechanics since the basic 
paper by Kay and Moses [I], who formulated the problem of constructing a plane 
stratified dielectric medium transparent to electromagnetic radiation in terms 
of a one-dimensional Schrodinger equation with a potential with preassigned 
bound-state spectrum that transmits without reflection continuum wave func- 
tions at all incident energies. From a mathematical point of view, the Kay- Moses 
method is equivalent to solving a non-linear Schrodinger equation whose poten- 
tial, V (x), is a quadratic function of a fixed number, n, of unknown bound-state 
wave functions it can also be considered as a kind of Hartree-Fock potential 
with n occupied states for a system of particles interacting through schematic 
contact interactions in one space dimension [3]. 

More recent approaches to reflectionless potentials in non-relativistic quan- 
tum mechanics make use, among others, of Darboux transformations 4 , super- 
symmetric hyerarchy derivations from the trivially transparent constant potential[5] 
and Casimir invariants of non-compact Lie groups [5], the latter method giving 
rise to large families of reflectionless potentials in implicit form, in addition to 
explicit analytical forms derived in previous approaches. 

In relativistic quantum mechanics, the Kay-Moses method has been applied 
to the one-dimensional Dirac equation with either scalar [7] [S] [5] or pseudoscalar 
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potentials [T^, since the presence of a vector component may break the trans- 
parency of the potential at all energies (see Ref. and Section 3 of the present 
work), with notable exceptions, one of which will be discussed in detail in Sec- 
tion 3. The relativistic extension of the Kay-Moses method is equivalent to the 
solution of an auxiliary non-linear Dirac equation. 

Reflectionless potentials play an interesting role in non-Hermitian theo- 
ries, too, such as quasi-Hermitian quantum mechanics [IT]. PT-symmetric quan- 
tum mechanics |1 2) . or pseudo-Hermitian quantum mecanics|14|.|15j. As is 
known, if a non-Hermitian potential V (x) is invariant under the product of par- 
ity V and time reversal T, so that VTV (x) {VTr^ = V* {-x) ^V{x), and the 
bound-state eigenfunctions of the Schrodinger Hamiltonian H — {x) 
are eigenfunctions of VT (exact VT symmetry) , the corresponding eigenvalues 
are real. As for the continuum of scattering states, it was proved in Ref. |16) 
for asymptotically vanishing potentials in the Schrodinger equation that if the 
asymptotic wave functions are eigenstates of VT (exact asymptotic VT sym- 
metry), the T^T-symmetric potential is reflectionless and unitarity is conserved. 
In Section 2 of the present work we extend the proof to the Dirac equation with 
potentials that admit non-zero constant limits at a; — >■ ±oo. 

Scattering from reflectionless potentials with exact asymptotic VT symme- 
try can thus be treated by the methods of standard quantum mechanics, with- 
out the need for an equivalent Hcrmitian formulation, which is neither exempt 
from technical difficulties, nor from ambiguity of interpretation: it has been 
shown that the equivalent Hermitian description of scattering from strongly lo- 
calized non-Hermitian potentials, a Dirac delta function with complex coupling 
strength in Ref. |17) and a T^T-symmetric combination of delta functions in Ref. 
[18] , implies strongly non-local metric operators and, consequently, an apparent 
breaking of causality due to incoming waves in the exit channel. This seems 
to be the price one has to pay in order to restore unitarity in the scattering 
process, although a new formulation of the problem|19| | 20 ) . based on the dis- 
cretization of the Schrodinger equation on an infinite one-dimensional lattice, 
has provided examples where the metric operator in the Hermitian equivalent 
formulation can be chosen as a diagonal matrix, called a quasi-local operator, 
which prevents the appearance of incoming waves in the exit channel, at the 
cost of a change of scale of the probability density on the left and the right of 
the scattering centre. 

In the present state of formulation of quasi-Hermitian theories, however, we 
share the opinion expressed in Ref.[l8j. that it makes sense to treat a non- 
Hermitian scattering potential as an effective one, accepting that it may well 
involve the loss of unitarity when attention is restricted to the system itself 
and not its environment, with which it can exchange probability flux (see also 
Refs. [21] [22] [23] ) . On the other hand, reflectionless potentials are a special 
class of 'PT'-symmetric potentials that conserve unitarity even in the standard 
formulation of quantum mechanics; therefore, we believe that they may deserve 
a study of their own, not only in the standard framework, i.e. with a trivial 
metric operator, adopted here, since it does not give rise to unphysical aspects 
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for an isolated system, but possibly also as a test of alternative approaches. 

The main scope of the present work is to investigate the behaviour of reflec- 
tionless potentials in relativistic quantum mechanics, under different conditions 
of Lorentz covariance, i.e. when they appear as vector, scalar or pseudoscalar 
components in the one-dimensional Dirac equation. Section 2 describes the 
general formalism, examples of scalar-plus-vector potentials are worked out in 
Section 3, pseudoscalar potentials in Section 4 and scalar potentials in Section 
5. Finally, Section 6 is dedicated to conclusions and perspectives. 



2 General formalism 

The time- independent Dirac equation in (1-1-1) dimensions with vector , scalar 
and pseudoscalar potentials, V, S and P, respectively, reads, in units h — c~l 

[axPx + /? (m + Six)) + laxPP {x) + V{x)] * [x) - E^^ (x) . (1) 

Here, ^' Ix) is a two-dimensional spinor and px = o-x and j3 are two 

ant i- commuting Hermitian traceless matrices with the property — (5"^ = 

^ ^ ; which can be identified with two Pauli matrices. 



1 

We assume for generality's sake that S, P and V can have non-zero limits 
at a; = ±oo: lim^-j-ioo S{x) = S±, and analogous notations for P and V. In the 

Dirac representation [24], ax = ctx = q ^ and (3 — (Tz — '^1 ) ' ^^"^ 

the asymptotic Dirac equation reads 

m + S±+V±-E -z(-|+P±) \f^Piix) 
^{-±+P±) -m-S± + V^-E J 

Let us search for a solution of the following form 



= . (2) 



ipi (x) = v4±e^'=±^ + S±e-*'^±"^ 
ip2 (x) = A±C±e*'=±^ + B±£i±e-*'=±^ 



(3) 



where A±, B± C± and D± are complex numbers. 

Direct substitution of formulae ([3]) in the asymptotic Dirac equation ([2]) 
yields 

fc± + iP± ^ -k±+ iP± 

^ m + S±+ E ~V± ' ^ m + S±+E-V± 
and the asymptotic momenta satisfy the relation 

kl^{E-V±)'' ~{m + S±)'' -Pi, (4) 

while A± and B± remain to be fixed on boundary conditions. 

Like in our previous works [^.126). which use the same representation of 
Dirac matrices, the parity operator, V, and the time reversal operator, T, are 

V^PoP = Po<7z , 

T — /C/3 — /Cctz , 
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where Pq changes x into ~x and JC performs complex conjugation. Their product 

VT - PolCcrl = PoJC (5) 

is thus the same as in non-relativistic quantum mechanics 16^. 

If the Dirac Hamihonian on the left-hand side of Eq. ([T|) is T'T-symmetric, 
S+ = S*_ , V+ = VI and P+ = -PI . VT symmetry of the potentials thus 
implies = fc^, which means that either k*_ = or fcl — —k+. Using 
formulae ([31), it is easy to show that only with the choice k*_ ~ fc+ the ratios of 
transmitted waves over incident waves remain finite even if the amplitudes of 
asymptotic wave functions may diverge at x = ±00 . This argument does not 
hold for the reflection coefficient: therefore, only when reflection is identically 
zero it makes sense to treat T'T-symmetric potentials that do not vanish asymp- 
totically. In turn, fc^. = fczp implies C± = and D^. = ZJzp. When all the 
potentials vanish asymptotically, the well-known expressions for free particles 
are recovered: fc^ = £'2 - m^, C± = -g^, D± = 

In general, A± and B± are linear combinations of the coefficients of asymp- 
totic expansions of two linearly independent solutions to Eq. ([Ij, ^^^^ (x) and 

Jm^ (a:) = a,± (A:±) ) e"^±- + (fc±) ) g-''^*^ = 1, 2) 

in the general asymptotic solution 

lim ^ (x) ^ a lim (x) -|- /3 lim (x) , 

x— >-±oo x— ^±00 ^±00 

or 

A± = aai± {k±) + /3a2± (fc±) , 
B± =Q6i±(fc±)-|-/?52±(/c±) . 

a and /3, in turn, can be fixed by boundary conditions. If ^I' (a;) is a pro- 
gressive wave, travelling from left to right, we must have, apart from a global 
normalization constant, not relevant in this context, 

lim,^_oo * (x) = J_ ^ e^'^-^ + Rl^r 

lim,^+oo * {x) = Tl^r ' 

where the transmission and reflection coefficients, Tl^r and Rl^r, have been 
introduced. 
Therefore 

= aai_ + /3a2- = 1 , 
B_ = abi^ + I3b2- = Rl^r , 

B+ = abi+ + pb2+ = , 
A+ = aai+ + Pa2+ = Tl^r , 
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whence 



_ fcl+a2+^ai+fc2+ 
L^R 6i_|_a2- — 01-62+ 
6i+b2-— fei-b2+ 



(6) 



bi + a2_-ai_b2+ 

In the same way, if 5* {x) is a regressive wave, travelhng from right to left 



hm,. 



^ {x) - Tr^l 



1 

^Rr- 



^ — ik- 



Thus 



whence 



A- = aai- + Pa2- = , 
B- = ah^ + pb2- - Tr^l , 
= aai+ + f3a2+ = Rr^l 
B+ = abi+ + Pb2+ = 1 . 

a2-hi- —01-62- 



61 + 02 



-01-62+ ' 



] D Q1 + O2- — Ql- 02+ 

R.^L 61 + 02-— 01-62 + 

Not surprisingly, the transmission and reflection coefficients 
the same as in the non-relativistic case|16|. 

The Wronskian of two solutions of the Dirac equation, 4''^^' {x) 

is defined as 



(7) 

and (III) are 

/ ^(1) {x) ^ 

V 4'' (^) ) 



and ^f^^) (a; 

W{x) E 





(^) \ 






(x) 




4'^ (x) 


4'' 



',(1) 



(2) 



',(1) 



(8) 



It is easy to check that — 0, i.e. W {x) — const., by expressing the 

derivatives of the spinor components as linear combinations of the components 
themselves, as dictated by Eq. ([T]). If the two solutions are linearly independent, 
W ^ 0, of course. 

Using definition ([8|) and asymptotic wave functions, 4'*^*''± = lim2;^±co (x), 
one easily obtains 

lim W (x) = W±^ (ai±62± - a2±fei±) {D± - C±) . (9) 

x—¥±oo 

Remembering expressions (jUlTl) of the transmission coefficients, formula ([9]) 
yields 

= (ai_62+ - a2-b,+)TR^L{D- - C_) , 
W+ - (01-62+ - a2-bi+) Tl^r{D+ - C+) . 

and — W+ is equivalent to 



Tr^l{D--C-)=Tl^r{D^ 



C+ 
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Tn^L D+-C+ D*_-C*_ ^ ' 

Here, v = 2 arg (Z3_ — C-) is a real phase. When, in particular, all potentials 
vanish asymptotically, C_ = — 1?_ are real numbers and the two transmission 
coefhcients are equal. 

The phase difference, of the two transmission coefhcients is different from 
zero when the imaginary components of the T'T-symmetric potentials do not 
vanish asymptotically and is present in non-relativistic quantum mechanics, 
too, as recently shown in Ref.|27) for a T'T-synimetric version of the hyperbolic 
Rosen-Morse potential. 

It is worthwhile to point out that the formalism just developed refers to local 
potentials. For non-local potentials it has been shown that the ratio of the two 
transmission coefficients is not 1, but a complex number of unit modulus, even if 
the imaginary potentials vanish asymptotically, both in non-relativistic [28^ and 
relativistic wave equations [25]. In this case the two reflection coefficients have 
the same phase, but different modulus and unitarity is broken. 

Let us now apply the VT operator ([5]) to the general asymptotic wave func- 
tions 

Jm^* (a;) = vl,± (x) =A±(^c±) + ^± ( D± ) ' 

or, more conveniently, to the following interpolating function, which coincides 
with the asymptotic wave functions at large \x\ 

(x) = ^(1 + sgn{x))-^+ {x) + ^(1 - sgn{x))^^ (x) . 

By definition ([S]) one gets 

VT'f^nt. {x) = i^x) = i(l - sgn{x))^l {-x) + i(l + sgn{x))^*_ {-x) . 

Imposing VT ^'mt. (x) ~ e'^'^'^int. {x), with ip a real phase, yields 

(-.t) = e'^P^^ (x) . 

Remembering the behaviour of k±, C± and D± under complex conjugation, 
we obtain the following constraints on A± and B± 

For a progressive wave (A_ = 1, = 0), this is equivalent to Tl^h = 
= e^^'^- and Rl^b. = B^ = 0, while, for a regressive wave ~ 0, 
B+ — 1), one obtains T^^l — B_ = e"'*^ and Rr^l = ^+ = 0. In other words, 
the potentials are refiectionless and conserve unitarity, since the transmission 
coefficients have unit modulus. 
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In the non-relativistic limit, \C±\, \D±\ ^ 1 and the lower components of 
Dirac spinors are negligible with respect to the higher ones. Non-vanishing 
potentials at x — ±cx) only affect asymptotic momenta k± and the preceding 
discussion and its conclusions remain valid, thus generalizing the case of short- 
range potentials treated in Ref.|16|. 

It is worthwhile to point out that potentials that behave asymptotically like 
T^T-symmetric step functions (P+ = _p* and so on) may admit asymptotic 
wave functions that are eigenstates of VT, unlike the step functions themselves, 
which are not reflectionless, because the asymptotic behaviour of wave functions 
is determined by the behaviour of the potentials in their whole domain. 

In the following sections, we specialize the general interaction of Eq. ([1]) 
to scalar-plus-vector, pseudoscalar and scalar potentials and, for each type of 
potential, work out some examples in detail. 



3 Scalar-plus-vector potentials 



In the present section, we specialize Eq. ([T|) to a scalar-plus-vector potential 
with the same x dependence: S (x) = csf (x), V (x) = cyf (x), with cs and cv 
real coupling constants. 

[a^Px + /3m + (cs/3 + cy) f {x)] ^ {x) = (x) , (12) 

We find it convenient to adopt the Dirac representation ax — Cx and l3 = az- 
The Dirac equation (IT^ satisfied by the spinor VI, (^) ^ Jig ^ is thus 
written explicitly in matrix form 

f m-E + {cv+ cs) f (x) -i± V ^1 = M 

V -m-E+ (cy - cs) f (x) )\^2(x) ) {O J 

(13) 

which reduces to a system of two first-order equations for the unknown spinor 
components -01 {x) and -02 (x) ■ In order to obtain analytic solutions, we limit 
ourselves to the particular cases cy = cs and cy — —cs, which correspond to 
spin symmetry and pseudo-spin symmetry in three space dimensions [29] . Let 
us consider the case cy = cs ~ c' first. It is easy to see that the two equations 
obtained from formula (jl3p reduce to the simple system 




(14) 

Here, the equation satisfied by ipi (x) is Schrodinger-like, with the same VT- 
symmetric form / (x) as the original Dirac equation and an energy- dependent 
potential strength s' (E) = 2c' {E + m) and tp2 (x) is obtained by deriving ipi (x) 
with respect to x. On the r.h.s. of the first equation, fc^ > for scattering states, 
while for bound states, fc^ < implies an imaginary value of k, corresponding to 
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poles of the transmission coefficient. In the limiting case fc = 0, both normaliz- 
able bound states and non-normahzable half-bound states [30], corresponding to 
transmission resonances, are possible, depending on the potentials under con- 
sideration. 

When cv = —cs — c", ■01 and '02 exchange their role, since '02 now satisfies 
a Schrodinger-like equation with the original / {x) and the energy-dependent 
strength s" (E) = 2c" {E ~ m), while 0i is proportional to the space derivative 

of '02. 

/ -ifj^2 (x) + 2c" {E - m) / {x) 2P2 (x) = {E^ - m^) {x) = PyJ2 (x) 

(15) 

Energy dependence of the coupling strengths in Eqs. (|14m5p may affect the 
reflection properties of a PT-symmetric potential. An example of this general 
behaviour is provided by the hyperbolic Scarf potential with integer coupling 
constants I and n 



P + n{n + l) 1 il{2n+l) sinhcc 
J^'''-' 2^ cosh^ X 2^ cosh^x ' ^ ' 



which is known to be reflectionless in the Schrodinger equation[T6] (note that the 
quoted reference uses units 2m — 1, as is common in non-relativistic quantum 
mechanics). When inserted in the Dirac equation, it gives rise to an equivalent 
Schrodinger-like equation (|14p where the potential maintains the same shape, 
but is no more reflectionless, because of the energy dependence of the coupling 
strengths. 

On the contrary, if / (x) exhibits an exact VT symmetry in the Schrodinger 
problem, it maintains it in the Dirac problem with the appropriate superposition 
of vector and scalar components, provided it is not connected with a particular 
value of the coupling strength, s' or s", which becomes a function of E. 

A notable example is provided by the T'T-symmetric potential 

/W = 7^-2. (17) 

(X + le) 

where e is an arbitrary real number, regularizing / at a; = 0, which is a well- 
known example of reflectionless potential in the Schrodinger case[T6|. Let us 
consider the case cy = cs = c' first, so that the equation ([T4)) satisfied by 'ipi 
reads 

- ^i^i (^) + ^^ + 7^ 1 (x) = {E' - m') V^i (x) ^ fc^^i (x) (18) 
dx (a; + ie) 

for scattering states (i?^ > m^). The above equation is Schrodinger-like and is 
quickly solved by introducing the complex variable z — k (x + ie) and factorizing 
'01 (z) = z^l'^Lp (z): in fact, the equation satisfied by Lp 



^{z)^Q (19) 



az'^ az 



z^ - 2c'(m + E) - - 
4 
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is a Bessel equation of index = 2c'{m + E) + j. Note that v is imaginary 
when 2c'(m + E) + j < 0, which can happen, for instance, for positive c' and 
large negative E, or viceversa. 

Two hnearly independent solutions of Eq. (jl9p with asymptotic behaviour 
appropriate to scattering states are the Hankel functions of first and second 
type, hI^^ (z) and hI"^^ (z), respectively 



lim|,|^,,i7W(^)^ (A)i/2exp[z(z-fi.-f)] , 
lim|,|^o,iji')(z)= (;^)'/'exp[-i(z-f;.-f)] , 



valid for SRi/ > —1/2, |argz| < tt, this latter condition being ensured by the 
non-zero imaginary part of z, i.e. 5z = ke. According to formulae (1141) . the 
corresponding linearly independent solutions of the Dirac equation are 



E-\-m dx 



where A = -gM — . In order to obtain the final form of the r.h.s., use has been 

made of the relation ■j^H^'''' (z) = hI'^\ (z) ~ ^U^v^ (2). The asymptotic be- 
haviour of the Dirac spinors is 

hm vl/(fe) [£) = lim ( ) , 

^-^±°o kKoo \^ (z) j 

or, more explicitly, using formulae (1201) 

Jm^ (x) ^(^^' ( 1 ) exp (zfcx - fee - ^\v - (21) 

and 



hm ^((2) (2;) ^ ^ ^ exp (^-i/cx-^fce-Hi^iz-f-i^) . (22) 



Formulae (I2m22|) are particular cases of the asymptotic formulae of Section 
2, whose constants now are 



1/2 

oicrt I — kf: — I 

2 4 

bi- = 61+ = 
0.2- = a2+ = 
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02- = 02+ = I ~ ) <3xp [ke + i—v + i- 
so that formulae immediately yield 

Tl^b. = Tr^l = 1 , Rl^r. = Rr.^l = . 

Of course, potential pT)) does not sustain bound states with fc ^ 0, because 
the transmission coefficients are independent of k and cannot have poles in k, 
or E. When fc = 0, a non-trivial solution can exist dX E ~ m: in this case, Eq. 
reduces to 

- (^) + T^T^^i = • (23) 

The solution to Eq. (|23|) can be searched for in the form of a power, (a; + ie)''^, 
thus leading to an algebraic equation for 7 

7^ — 7 — 4c'm — , 

whose solutions are 



1 ± (1 + 16c'm) 
71,2 = 



1/2 



2 

Depending on whether l + lGc'm ^ 0, the two roots are either real or complex 
conjugate: in both cases, the general solution to Eq. can be put in the form 

ipi (x) — ai{x + iey^ + a2{x + iep^ , 

where Ui (i = 1,2) are to be fixed on boundary conditions. It is easy to un- 
derstand that a normalizable solution, i.e. a bound state, can exist only when 
1 -|- 16c' m > 1, or c' > 0, by choosing ai = 0. In this case, the solution for 
ipi (x) reads 

ipi{x) =a2{x + ie)~ . (24) 



Here, /? — \/l + 16c' m > 1 and 02 can be determined by normalization of 
the complete Dirac spinor 

-t-00 

/ dx ( r^ (x) ^±^1 ix) ) ( __^^(^| {x))=^- (25) 



Both integrals in formula (|25p can be computed analytically in terms of 
asymptotic expansions of the hypergeometric function, F [A, B, C; z), since they 
can be reduced to the integral representation 

J {l + y^y"dy^yF(^^,a,^-,^i/^ + const, (^^a > ^ 

yielding 

'•+00 r(rj-i) 
(1 + y) dy ~ y/n 



r(a) 
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where T (a) is the Euler gamma function. 

The final result for the normalization constant is 



ml = 



1 

4m^ 



2 



For the sake of completeness, we mention the case of the double root, 7 = 
1/2, of the characteristic equation, occurring when c' = —j^- the general 
solution to Eq. for this case can be put in the form 



-tpi (x) — ai {x + ie)^^^ + a2 (a; + ieY'^ In (x + ie) 

and is not normalizable. 

The case cy = —cs = c" can be treated in a similar way starting from Eqs. 
(fT5|) : the solution of the Schrodinger-like equation for ip2 (x) is obtained in the 
same way as before and the bound state with fc = now appears at E = —m. 



1/2, 



4 Pseudoscalar potentials 

The Dirac equation with a pseudoscalar potential, P (x) = cpf {x), and cp a 
real coupling constant, reads 

[a^Px + 13m + icpa^Pf (x)] * (x) = (x) . (26) 
In the Dirac representation, iUxP — i<Jx<^z = "'y ~ ( ^ ) ■ invari- 



X * ^ 

ance of the interaction term implies /* {—x) = —f{x). After expressing Eq. 
([26]) as a system of coupled equations in the two components, V'l (x) and 1^2 (x), 
of the Dirac spinor, ^ (x), 



m~E ~i±~^P (x) \ ( ^Ax) \ _ ( 



i^+iP{x) -m-E )\i})2{x)) \Q 



(27) 



it is almost immediate to derive the two decoupled Schrodinger-like equations 
satisfied by ■0i and -02 

(-^^ + U,{x)^ ^0 (x) = (x) = e^, (x) , (j = 1, 2) (28) 

where U^{x) ^ ^{PHx) + ±P{x)), U2 (x) ^ ^{P' (x) - ±P{x)). As 
already shown in Ref.|31|. the two T^T-symmetric Hamiltonians 



i^.^-^^ + C/.(.U.-l,2) (29) 



constitute the Bose sector of a non-Hermitian representation of an s/(l|l) superalgebra 
The corresponding super-Hamiltonian is 
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while the differential operators 



forming the two partner Hamiltonians Hi — LM and H2 = ML, give also rise 
to the "supercharges" 

which form the Fermi sector, since H, Qi and Q2 are closed under the following 
set of commutation and anticommutation relations 

{Qi,Q2}+=n, [Qi,n]_ = [Q2,n]_ = o . 

It is worthwhile to point out that, differently from the Hermitian case, Q2 
is not the Hermitian adjoint of Qi, but the following relations hold 

M = {VT)L {vry^ = vL^v-^ , Q2 = vq\v~^ . 

One thus speaks, in this case, of a 7^-pseudo-supersymmetry . Note that 
P(x) plays the role of a superpotential. The fact that the supercharges ([5^ 
are not the Hermitian conjugates of each other gives rise to a richer structure 
of supersymmetric systems compared with Hermitian theories [33] . 

As is known, the discrete spectra and the scattering properties of the partner 
Hamiltonians are connected by supersymmetry, so that it is sufficient to compute 
bound states and scattering states of one partner only. 

In particular, it is not difficult to prove the relations connecting the trans- 
mission and reflection coefficients of the two partners, in terms of the asymp- 
totic limits of the superpotential, P± = \iinx^±Qo P {x) and the corresponding 
asymptotic momenta, k"^ = E"^ — m? — : let T^}^j^ (RlI^r) be the transmis- 
sion (reflection) coefficient of a progressive wave in presence of potential ?7*-^^ 
(j = 1, 2) : we easily obtain 

_ -»fc_+P_ n(2) 

Similar relations connect the transmission and reflection coefficients of a 
regressive wave 

p(l) _ zk++P+ p(2) 

^(1) _ -ifc__KP_ ^(2) y.'^^) 

^R^L — -ik++P+^R^L 

Derivation of formulae p3ll34p is given in Appendix , for the sake of com- 
pleteness. 

It is obvious that if one of the two partners is reflectionless, so is the other 
and, consequently, the superpotential in the Dirac equation. 
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In Ref.[Tn], the Kay-Moses method of constructing reflectionless potentials 
was extended to real symmetric pseudoscalar potentials in the Dirac equation. 
The examples worked out in the quoted reference can be made PT-symmetric 
by applying an imaginary shift to the space coordinate, x ^ x + ie. Thus, the 
following superpotential, with cp ~ 1 for simplicity's sake 



/ (x) — tanh (x + ie) + 



1 



tanh {x + ie) — Acoth (A (x + ie)) 
and A > 1, generates the following supersymmetry partners 



(35) 



2m 



and 



A2 



2 (A^ - 1) (A2 cosh^ {x + ie) + sinh^ (A {x + ie))) 
(A cosh (x + ie) cosh (A {x + ie)) — sinh (x + ie) sinh (A {x + ie)))^ 



U2{x) 



1 
2m 



(36) 



cosh {x + ie) , 

A part from the constant term }? j (2m), which enters in the definition of the 
asymptotic momentum, [x) is a reflectionless Poschl- Teller potential[T2] and 
V\{x) is necessarily reflectionless, too. It is worthwhile to point out that our 
definitions of IJ\ and U2 are exchanged with respect to Ref.[10], but in agreement 
with Ref.[3T]. Transmission and reflection coefficients for potential p6|) can be 
immediately written down from the corresponding formulae of the more general 
hyperbolic Scarf potential obtained in Rcf.[16 , after observing that /± — =fA 

and the asymptotic momenta are k± - 



(2) 



R 



£:2 



- ^R^L - ^ ' 
^R.^L l + ik 



l2 - A2 = fc 



(2) 



T 



(2) 

R.^L 



(37) 



1, equivalent to 



A real fc is a necessary condition for 

exact VT symmetry of the asymptotic wave functions. 

Formulae (1331134^ immediately yield the corresponding coefficients for poten- 
tial Ui 



R 



(1) 



(1) 

L-S.7? 



R 



(1) 



, 



(1) 



R-^L 

\-ik 1 



'L-j-i?, ^R-^L \+ik l+ik ■ 

As far as bound states are concerned, it is well known [5^ that U2 {x) admits 
only one bound state with eigenvalue e = ( or i?^ = ^2 _|_ ^2 _ gj^^ ^]^q 
corresponding wave function is 

V'2 {x) 



N 



cosh [x + ie) 



(38) 



with the constant N to be determined from normalization of the complete Dirac 
spinor. 
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Since i?2'02 (a;) = MLip2{x) = ^;^V'2 (a;), we have that LH2^J2{x) = 
LMLip2 (x) = i?iiV2 (x) = ^^Lip2 {x). Therefore Vi {x) = coLtp2 (x), with 
Co a normahzation constant, is eigenfunction of Hi with eigenvalue 2m ^^'^ 
corrresponds to the first component of the Dirac spinor. From the definition of 
the differential operator L and from the first equation (|27p we get cq = ^^^^ 
and 

iN \ A2 - 1 



— E J sinh (a; + ie) — A cosh (x + ie) coth (A (a; + ie)) 

+00 +00 



= 1. 



with / ^-t (a;) * (a;) / f iV'i (x)| V |V'2 (x)n dx 

— CO —CO 

V'l (x) from formula p9p has a node at x = — ie and is not the ground 
state of Hi . Conversely, if a non-trivial normalizable solution of the equation 
M-ipi (x) — exists, then (x) = LM^pi (x) — 0, and tpi (x) is eigenstate of 

Hi with eigenvalue e = 0. In this case, -0]^ (x) cannot be written as Cf)L'ip2 (x), 
with 11)2 (x) a non-trivial normalizable function, since, otherwise, we would have 
ML'02 {x) = i?2'02 (^) = and t/ij (x) would be an eigenstate of H2 with eigen- 
value £ = 0, which is impossible, because H2'4'2 ix) = admits only the trivial 
solution ■02 (x) = 0. 

The equation 

{x)^^l^-j- + f{x)^liii{x)=Q, (40) 

is satisfied by 

ipi (x) = exp I / / (x') dx' 



Note that the condition limj;^±oo / (x) = /± = T^, with A > 1, yields 
lima;^±oo -01 (x) = 0. 
In our case 

j f (x) dx — In (coshx) — In (2A cosh x cosh (Ax) — 2 sinh x sinh (Ax)) -|- InTV , 
so that 

- N 

w,; (x) — (41) 

^ ^ Acosh(A(x + ie))-tanh(x-f ie)sinh(A(x + je)) ' ^ ^ 

where A'' is to be determined from normalization. The second component of 
the Dirac spinor, 02 (x), is solution of the equation H21P2 (x) = ML'iJj2 (x) = 



14 



^ Lip2 (^) = 0) which admits only the trivial solution, {^) =0- iV is thus 

+00 +00 2 

determined from the condition J {x) [x) dx ~ J \tpi(x)\ dx — 1. 

—00 —00 

— 1 

In this case, the ground state of H2 has the same energy, e = > 0, as 

the first excited state of Hi, whose ground state has e = 0. This is an example 
of exact supersymmetry. 

A second example of reflectionless pseudoscalar potential, whose bound states 
were already studied in Ref.|31j. is 

f (x) = ntanhx + i — — , (42) 
cosh a; 

with integer constants n and I. In fact, the two supersymmctric partners from 
formula (|29p are 

Ui {x) = ^(n^- + il (2n - 1) -^i^ 

^ \ ' 2m \ cosh^ X ^ ' cosh'^ x ^ 

U2 (x) = ^ ( n2 - !i(n+10±£ + (2„ + 1) sinh. 

\ > 2m y cosh^ x ^ ' cosh^ £ 

2 

which, a part from the constant term are reflectionless potentials of hyper- 
bolic Scarf type (|16p . In this case, f± = ±n and k± = \JE'^ — m? — = k. 
Here again, the transmission coefficients are given by Ref. [16] 

^(2) ^^(2) _ ( -. y»+i {n~ik)...{l- ik) {I - \ -ik) ...{\- ik) 
for n > 1 and 

^(1) ^^(1) _ , . .n+i-i {n-l-ik)...{l-ik){l-\-ik) ...{\~ ik) 

^L^R ^R^L I ^) {n-l + tk)...{l + lk){l-l+ik)...{^ + ik) 

for n > 2. 

Unitarity and asymptotic VT symmetry are conserved if k is real. 

As for bound states, U2 (x) admits n of them, all with real energies, and 
Ui (x) has n — 1 bound states at the same energies of those of U2 (x), excepted 
the ground state of the latter. Here again, the pseudo-supersymmetry is exact. 



5 Scalar potentials 

When only a scalar potential, S (x) = S* {—x), is present, Eq. ([ij simplifies to 

[axPx + I3{m + S {x))] * (x) = (x) . (43) 

In spite of its apparent greater simplicity, however, Eq. (|43l) is more difficult 
to solve than Eq. (jl2p . including a vector potential with the same x dependence 
of the scalar potential and equal, or opposite coupling strength, if one adopts 
the Dirac representation — (Tx, (3 — cTz, because the second order equation 
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satisfied by the first component, ij^i (x), of the Dirac spinor, {x), now contains 
also a first order derivative, -^ipi (x), whose x-dependent coefficient is negligible 
only at a; = ±00, provided S (x) admits constant limits, lim2;_>.-|-oo S (x) = S±, as 
assumed in Section (2) in deriving the general form of asymptotic wave functions 
for arbitrary combinations of scalar, vector and pseudoscalar potentials. 

If we are interested in exact wave functions, including those corresponding to 
bound states, it is more convenient to adopt a different representation of Dirac 
matrices, ax = cry, (3 = ctx, like in Ref.(31). The main drawback of this choice 
is that the kinetic term of the Dirac Hamiltonian, axPx = cryPx, is not VT- 
symmetric, but the remedy is simple, since the two representations are unitarily 
equivalent: the unitary transformation 

U = e^^^'e^^'^^ = l[l + ^{ax+ay+az)]=^-i^\+^^ |^.) (44) 

changes the and (3 matrices of the representation of Ref. [JT] into those of the 
Dirac representation, since 

UdyU^^ = ax , 
UuxU-^ = . 

It is also of interest to determine the operator that corresponds in the 
present representation to to the VT operator ([5]) of the Dirac representation: 
let us notice that VT anticommutes with iax, since T is antilinear, and that 
(VTia-x)^ = I2, the 2x2 identity matrix. It is not difficult to check that 

UVTiaxU-^ = -iUaxU^VT = VT . 

Therefore, VTifJx is obtained from VT by means of the similarity transfor- 
mation that connects matrices in the present representation with those in the 
Dirac representation. Moreover, once the Dirac equation (|43)) has been solved 
in the new representation, and the spinor ^ [x) is known, the corresponding 
solution in the T^T-symmetric Dirac representation will be 

-^oix) = C/*(x) . (45) 

As a consequence, if d [x) is an eigenstate of VT, 4* (x) is an eigenstate of 
VTiax- 

In the new representation the two equations satisfied by the components, 

ipi {x) and ■(/'2 [x), of 5* (x) 



-E -d_ + m + S{x) 

l+^ + S{x) -E 

are easily decoupled as 



i>2 (X) I 



- 2^^^! W + 2^ ((™ + ^ - - 1-xS (x)) ^1 (x) = ^V-i (x) ^ e^i {x) 

(x) = -i(£ + m + 5(x))^i (x) 



(46) 
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where the equation satisfied by ipi (x) is Schrodinger-hke, with an effective po- 
tential 

Ui ix) ^±-(^im + S ix)f ^rn^^^S (x)^ . (47) 

If one derives instead the equation satisfied by ip2 (x) , the resuh is a Schrodinger- 
hke equation with the same effective energy e = ^ ""^ and an effective potential 



U2{x)^^[{m + S{x)y 
Zm 



(48) 



We are thus led again to a pseudo-supersynimetry 31 , like in the case of a 
pseudoscalar potential. Note that S (x) + m = m (x), the effective Dirac mass, 
now plays the role of a superpotential. The supercharges now are 



L = 



'2m 



dx 



+ m + S {x) 



M = 



1 



/2rri \dx 



Six] 



and the partner Hamiltonians 



In this representation, the asymptotic Dirac equation 



_d_ 

dx 



-E 

- m + S± 



_d_ 

dx 



- m + S± 
-E 



■fAii [x] 

■>p2± (x) 



where, as before, lim^^±oo (x) = S±, limx^±oo ipk (x) 
solved in the form 



ipk± (x), is easily 



*± (x) 



A. 



tpi± (x) 
i'2± (x) 
1 

C± 



Jk±x 



B± 



1 



-ik-^x 



where k± = I E — (m + S± ) ) are the complex asymptotic momenta and 



E 
E 



while A± and E± are to be fixed on boundary conditions. Owing to the fact that 
^± {x) and '^D± (x) are connected by Eqs. (|44II45I) . the asymptotic expansion 
coefficients A± and B± are related by the inverse of formula (|45|). i.e. ^ (x) = 
t7^5'£) (x), to the corresponding ones in the Dirac representations, A± and B±, 
given in formulae ^ with V± = P± = 0, in the following way 



A±^-i2L m+S±+E-ik± 
1 



m+S±+E 
m+S-|-+_E+ifc-|- 



(49) 
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It is worthwhile to stress once again that in the present representation of the 
Dirac equation 'i>± {x) is not eigenstate of VT- 

Before discussing specific examples, it is worthwhile to recall that a reflec- 
tionless potential can be obtained as a supersymmetry partner of the constant 
potential: in fact, putting Ui (x) = c in formula (|47p . one obtains a Riccati 
equation for the superpotential S (x) + m, which can be solved by separation of 
variables, giving rise to different solutions in connection with the sign of the con- 
stant m+2c. The superpotential is a trigonometric function of x when m+2c < 
(or c < — y) ^^'^ 9, hyperbohc function when m + 2c > 0; consequently, U2 {x) 
is a trigonometric Poschl- Teller potential in the former case and a hyperbolic 
Poschl- Teller potential in the latter. In the limiting case m + 2c = 0, one obtains 
S {x) + ni = —5;^, where d is an arbitrary constant, and U2 (x) = ^ (^^.^cif ~ T' 
When d — ie, with real e, U2 {x) is qualitatively similar, apart from the additive 
constant, to the potential studied in Section 3. Similar considerations could 
be made in the case of the supersymmetry involving pseudoscalar potentials 
presented in Section 4. 

A simple example of reflectionless scalar potential is the PT-symmetrized 
form of the real potential with one bound state derived in Ref. 

= F Si 1 ^ (50) 

_4 \ _ 

Eb cosh {kb (a; + ie) — As) cosh (k^ (x + it) + As) 

Here, Eb = and kb ~ \J rn? — E^ — are energy and momen- 

tum of the bound state, expressed as functions of the coupling strength C5 in the 
auxiliary non-linear Dirac equation discussed in the above mentioned reference 
and \b — |arccosh^-^^ . The two partner potentials from formulae (|47ll48p 
can be written in the compact form 

^ y ' m cosh^ (KB^2;+ic)-A_B) ' (51) 

^ ^ ' m cosh'-^ (kb (x+ie)+XB) 

While S (x) is T'T-symmetric, the partner potentials ([^Tjl are not, as ex- 
pected from Eqs. (|47ll48p . The two Schrodinger-like equations with potentials 

m 

1 d^ 

- 7^-^^■^^{x) + U^{x)^|^^{x)^ei}^{x) (i-1,2) (52) 
are satisfied by the bound-state wave functions 



Ni 



V'l (-^) cosh(K.B (a;+ie) — Ab ) 
'^2 (^) — cosh{KB{xXie)+XB) 



(53) 



respectively, where iVi and N2 are normalization constants. Note that the two 
partners have the same discrete spectrum, i. e. one bound state with real energy 
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s — the pseudosupersymmetry is thus spontaneously broken [5]. The two 

normahzation constants are related by the second Eq. (j46l) : after replacing in it 
formula (ISUI for S {x) and the first formula (1551) for ipi (x) , some straightforward 

■01 (x) 



algebra gives N2 = Ni. The corresponding Dirac spinor, (x) — 1 ^ 

has real energy Eb = y^m'^ — k\ and is normalized to one. Hence, 

+00 +00 +00 

■^^x)-^{x)dx= I \^x{x)\^dx+ [ \ij2{x)fdx 



+00 

4 ,2 /• dy 4e|7Vi|' 



|7Vi 



KB J cosh {2y) + COS {2k Be) sin (kbe) cos (kbe) 

—00 

2 sin (K^e) cos (Kse) 



4e 

Eqs. ([52]) with e > are satisfied by the scattering wave functions 



(54) 



respectively, with fc = v2me, if boundary conditions for incident progressive 
waves (L — ?> i?) are imposed. Since limj:_>.±oo tanh (Kf, (x + ie) ± As) = ±1, 
'(/'z (a;) does not contain a reflected component, so that Ri^l^n — 0. It is also 
immediate to determine the transmission coefficient 

_ lim^r^+oo ipi jx) _ ik - Kb 
hm:r^_oo ■04 {x) ik + Kb 



Since Tj^\j^ = Tj^}^j^, formula ([55)) yields the transmission coefficient T^^n 
of the Dirac spinor, which has unit modulus, as expected. The replacement k — 
—k changes progressive waves ([5^ into regressive waves, whose transmission 
coefficient now is 



^ lim^^_oo ip, (x) ^ -ik + Kb ^ ^(t) 

lim^^+00 V* (a;) -ik-Kb ^^-^ ^ ^ 



with = R^l^R — 0- The equality of the transmission coefficients could 



have been proved also by means of formula (fTO|) . since, in the present case, 
C_ = C+ = and D^=D+ = ^IM™. 



6 Comments and outlook 



In the present work we have extended the analysis of reflectionless PT-symmetric 
potentials in non-relativistic quantum mechanics presented in Ref.^lGj to rela- 
tivistic quantum mechanics with different forms of potentials in the Dirac equa- 
tion: scalar, pseudoscalar , or a mixture of scalar and vector potentials. We 
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have examined the connection between reflectionlessness and exact VT symme- 
try of bound-state wave functions and asymptotic wave functions even in the 
case the potentials have non-zero hmits at a; ^ ±00, thus removing a constraint 
imposed in Ref.|16j. Along this line, one could study further reflectionless po- 
tentials that diverge at a; — > ±(X) , such as the class of real symmetric potentials 
V {x) — —x^^^"^ [k = 1,2,...) discussed in Ref . |34) . A main drawback is that 
they are not exactly solvable in general and one has to resort to some approxi- 
mate method, such as the WKB method of Ref. j34j . 

Reflectionless potentials are expected to play a peculiar role also in VT- 
symmetric quantum mechanics in higher dimensions, very little explored up to 
the present time and only in non-relativistic problems. In three dimensions, for 
instance, only non-central potentials can exhibit non-trivial VT symmetry: in 
polar coordinates r, 0, 0, they satisfy the relation[35] 

V{y) = V (r, e, 0) = T/* (r, vr - 61, + tt) = V* (-r) . 

General characteristics of bound states are discussed in Ref. [3 5) in case of 
exact VT symmetry and in Ref. |i36) when the symmetry is spontaneously broken. 

Scattering states of PT-symmetric potentials in higher dimensions have not 
been discussed so far, but a complete analytic description should be possible if 
these non-central transparent potentials can be related to the euclidean group 
in n dimensions, E{n), for n > 3, since it is the maximal symmetry group of 
the transparent null potential in n dimensions. For instance, Ref. [371 solves two 
classes of real transparent potentials in three dimensions that admit E (4) as 
the potential group, in the sense that the corresponding Hamiltonians depend 
on the restriction of the quadratic Casimir operator to subspaces appearing 
in two subgroup reduction chains of E{4). The possible group-theoretical 
treatment of complex potentials would require non-standard realizations of the 
non-compact Lie group involved (non- unitary representations), along the lines 
of Ref. [38j . Extension of this study to relativistic quantum mechanics would be 
possible by relating to the Casimir operator of a non-compact group the mass 
invariant operator of the system within the framework of a Bakamjian- Thomas 
realization of the generators of the Poincare group, under which one has to 
ensure the invariance of the scattering matrix 

Finally, it is worthwhile to mention that interest in the first quantized ver- 
sion of the Dirac equation has been renewed in recent years by quantum sim- 
ulations of the dynamics of Dirac fermions with controllable laboratory sys- 
tems underlying the same mathematical models: single trapped ions are par- 
ticularly suited to this purpose and experiments have been designed to simu- 
late the one-dimensional trembling motion (Zitterbewegung[40] ■ |41) ) of a free 
Dirac fermion[^.|43j. the overcriticality effects of a vector potential well in one 
dimension |42|. the relativistic Landau levels generated in a constant homoge- 
neous magnetic field [44] . the pseudoscalar Dirac oscillator in two dimensions |45|. 
The possibility of simulating in quantum optical systems the one-dimensional 
Dirac dynamics in presence of non-Hermitian potentials, such as the T'T-symmetric 
ones considered in the present work, would be of mutual benefit to both fields 
of research. 
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A Relations between transmission and reflection 
coefficients of supersymmetry partners in the 
pseudoscalar case 

Let us consider a scattering solution of one of the two partner Hamiltonians, e. 
g. Hi, with energy e = ^ ^7 " > 

(x) = £0(1) (x) 

and a solution, 4'^^\ of Hamiltonian H2 with the same energy e 

7720(2) (x) = £0(2) (x) . 

Since Hi — LM and H2 — ML, where L and M are the differential oper- 
ators defined in formulae ([5T|) . we immediately see that £L0(2) = Li/20(2) = 
LML0(^) = iJiL0(^), which means that L0(^) must be proportional to 0(i), or 

0(1) (x) = CL0(2) {x) , 

where C is a constant to be determined. Let us assume now that 0(2) (x) is a 
progressive wave, with asymptotic behaviour 

lim,^_oo 0(2) (x) = e*'^-- + E'lli^e-''^-^ , 
lim,_,+oo 0(2) (x) = rl2_^fle^'=+^ , 



where fc± = ^2to (e — C/2 (±00)) = y — m"^ — P| are the asymptotic mo- 
menta. Thus, we have 

(1) (x) = Clim,^_ooi0(2) (a:) = ^ (£ + P_) (e^*=-- + 4'-^fle"''-") 



c 

/2m 



and 

lim,^+oo 0(1) (x) = Clim,^+eo £0(2) {x) = + P+) rf2^^e^'=+- 

On the other hand, we know that 

0(1) (x) = e'^-'' + Rl^rG"''"-'' , 
hm,^+oo 0(1) (x) = rl'^^e*^-+- . 
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Hence we obtain formulae (1551) of the text 

C _ 1 



^2m ik^+P- 
_ -»fc_+P_ p(2) 

rj,{l) _ ik++P+ rp{2) 

^L^R — ik^+P-^ L^R 



One proceeds in a similar way for regressive waves, 0^^-' (x) = C'Lc? 
with asymptotic behaviour 



lim,^_o, (x) = r^'2,ie-^^-- , 
lim,_,+oo (a;) - + iJ^l^ie^'^^- 



a = 1,2) 



with the result 



nU; _ TDK 
^^R^L — -ik+ + P+^^R^L , 
^(1) _ -ik_+I\_rp{2) 
^R^L — -ik+ + P+ R~^L 

corrresponding to formulae (p4l) of the text. 
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